Multi-scale Correlation Functions in Strong Turbulence 
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Under the framework of V. Yakhot [Phys.Rev.E, 57, 1737 (1998)] modelling of intermittent structure 
functions in fully developed turbulence and based on the experimentally supported Markovian nature 
of turbulence cascades [R.Friedrich and J.Peinke, Phys. Rev. Lett, 78, 863 (1997)], we calculate the 
multiscaling correlation functions. Fusion rules [V.S.L'vov and I. Procaccia, Phys. Rev. Lett. 76, 
2898 (1996)] which are experimentally tested [R. Benzi, L. Biferale, and F. Toschi, Phys. Rev. Lett. 
80, 3244 (1998)] to be compatible with almost uncorrelated multiplicative process can analyticaly 
be supported by direct calculations. 

PACS numbers, 47.27.Ak, 47.27.Gs, 47.27.Jv and 47.40. Ki. 



One of the most important issues in stationary turbu- 
lence is the intermittent behaviour of the velocity fluctu- 
ations in the inertial range. Understanding the statistical 
properties of intermittency is one of the most challenging 
open problems in three dimensional fully developed tur- 
bulence. The structures that arise in a random flow of 
stationary turbulence resembel themselves as high peaks 
at random places and random times. The intervals be- 
tween them are characterised by a low intensity and a 
large size. Rare peaks are the hallmarks of PDF's non- 
guassian tails. These strongly non-Gausssian activities 
are statistically scale-invariant processes responsible for 
energy transfer. Intermittency in the inertial range is 
usually analysed by means of the statistical properties of 
velocity differences, 5ru{x) = u{x -f r) — u{x) The 
overwhelming experimental and theoretical works have 
been brought forward for characterisation of structure 
functions; i.e. Sp — {{Sru{x))P) . A wide agreement exists 
on the fact that Sp{r) exhibits a scaling behaviour in the 
limit of high Reynolds number, that is Sp{r) ~ (x)'''' 
L » r >> rjk , where L is the scale of energy injection, 
rjk — {^Y^^ is the dissipative scale, e is the mean energy 
dissipation range and v is the kinematic viscosity. Rare 
peaks in the random flows are signaled with a nonlinear 
form of C,{p). In other words the velocity increments are 
multifractal and C(p)'s do not follow the celebrating K41 
theory, C,{p) — p/3. Recently it has been proposed 
that it would be more natural to look at single time corre- 
lations among velocity increment fluctuations at different 
scales: 

Tnix\ri,r2, ...,rn) = {5r^u{x)5r^u{x)...5r„u{x)) (1) 

where all the scales are lying in the inertial range, i.e. 
rj « ri « L For simplicity we are confining the dis- 
cussion to longitudinal velocity increments. Fusion rules 
that describe the asymptotic properties of n-point 
correlation functions when some of the coordinates tend 
toward one other are derived from two fundamental as- 
sumptions which are of paramount importance for the 



description of non-perturbative aspects of the analytic 
theory of stationary turbulence. The fusion rules are 
tested experimentaly and a good agreement between ex- 
periment and theory have been observed Q . li p < n 
pairs of coordinates of velocity differences merge, with 
the typical separations ri r for i < p and the remain- 
ing separation at the order of R such that r « R « L 
, the fused multiscale correlation is deflned as 

Tp+g{r, R) = {{u{x + r) - u{x))P{u{x + R) - u{x)y) 

= {iSrUix)n6Ruix)r) (2) 

it has been deduced that 



Tp+,{r,R) ^ Spir)Sp+,iR)/SpiR) 



(3) 



On the other hand recently multiscale correlation func- 
tions in high Reynolds number experimental turbulence, 
numerical simulations and synthetic signal was investi- 
gated by Benzi and coworkers |^ and it has been found 
that whenever the simple scaling ansatz based on the un- 
correlated multiplicative processes Q is not prevented by 
symmetry arguments, the multiscale correlations are in 
good agreement with the fusion rules prediction even if 
strong corrections due to subleading terms are seen for 
small-scale seperation r/i? ~ 0(1) . All the findings has 
led to the following conclusion that multiscale correla- 
tion functions measured in turbulence are fully consis- 
tent with a multiplicative almost uncorrelated random 
process for evolution of velocity increments in scale . Al- 
though successful interpretation of the fusion rules can be 
realised by considering a multiplicative random process 
for the evolution of velocity increments in length scale, 
but it is at most a phenomenological modelling and it is 
not based on a flrst principle calculations. Other front of 
experimental investigations over the bahaviour of condi- 
tional probability densities of velocity increments in scale 
have shown that Markovian nature of velocity increments 
in terms of length scale and in the inertial range would 
support the experimental data |7j. In fact the necessary 
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condition of " Markovianity" for velocity increments have 
been tested experimentaly from which phenomenological 
senarios for modelling the intermittency have been de- 
veloped 0. Afterwards the aforementioned ideas were 
supported by invoking to the theoretical ideas inspired 
by Polyakov Q and Yakhot [|j based on OPE and gen- 
eral invariances of the Navier-Stokes equation |^ . In this 
paper we show that the interpretation given by Benzi 
and coworkers Q dictating an uncorrelated multiplica- 
tive process for velocity increments in the inertial range is 
not just a mere phenomenological interpretation and can 
be supported by Yakhot modelling. Even we will candi- 
date the underlying dynamical process in scale which in- 
corporates the fusion rules of multiscale correlation func- 
tions. The calculations are consistent with the picture 
of almost uncorrelated random multiplicative process at 
least in the Fokker-Planck description of turbulence cas- 
cades. Furthermore we are able to connect the fusion 
rules to the Markovian nature by a simple operator for- 
malism even by perserving all the terms in the Kramers- 
Moyal's evolution operator of velocity increments. 
Let us start with the Navier-Stokes equation: 



vt + (v • V)v = j/V^v 



P 



+ f(x,i), V-v = (4) 



for the Eulerian velocity v(x, t) and the pressure p with 
viscosity v, in N~dimensions. The force f(x, i) is the ex- 
ternal stirring force, which injects energy into the system 
on a length scale L. More specifically one can take, for 
instance a Gaussian distributed random force, which is 
identified by its two moments: 

(/^(x,t)/,(x =fc(0),5(<-<>^,(x-x) (5) 

and (/^(x, t)} — 0, where ^,,v — xi,X2, - ■ ■ ,xn. The 
correlation function k^^(r) is normalized to unity at the 
origin and decays rapidly enough where r becomes larger 
or equal to integral scale L, that is k^^{rij) = ^(0)[1 — 



c _ (ri,)^(r,j), 
-Xil 



] with fc(0) and L = I, where r.; 
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Recently Yakhot generelize the Polyakov's approach 
on the Burgers turbulence for strong turbulence. He 
used the OPE approach to closed the equation for the ve- 
locity increment PDF and show that the PDF of longtu- 
dinal structure function S*, — {{u{x + r)~u{x))'^) ~ {U'^), 
where u{x) is the a;-component of the three-dimensional 
velocity field and r is the displacement in the direction of 
the a;-axis satisfy the foUwing equation for homogeneous 
and isotropic turbulence in the limit r — > as, 



_d_^dP _^ dP_ _A_d_ 
dU dr ^ dr r dU 



UP 



L at/2 



UP (6) 



where A 



B = -So > and for the Navier- 
Stokes turbulence it has been shown that i? ~ 20 can be 
derived by a self-consistent calculations j|] . Last term in 



the right hand side is responsible for the breakdown of 
Galilien invariance in the limited Polyakov sense, which 
means that the single-point Urms induced by random 
forcing enter the resulting expression for velocity incre- 
ment PDF's. 

Now one can show that the probablity density and as 
a result the conditional probablity density of velocity dif- 
ference satisfy the Kramers-Moyal's evolution equation 

Eil: 



dU' 



-{D(^\r,U)P) 



(7) 



where D^'^\r,U) = ^^[/" -I- /3„C/"~^ We have 
found that the coefficients Q!„ and /3„ depend on A 
and B, u^ms a-nd inertial length scale L which are 



given by a„ = (-1)'^ 



and — 



(B + l)(B+2)(S+3)---(S+n) 

(-1)" "T° (B+2){B+3)...(B+n) ^liere /3i = by homogene- 
ity Q. The coefficients D'^''\r,U) are the small scale 
limit of the conditional moments They fully char- 

acterise the statsitics of eddy distribution in the inertial 
range which is defined as: 

hm,,^,, ^J{Ui - U2rP{Ui,n\U2,r2)dUi 

It is noted that P{Ui,ri\U2,r2) statisfies the equation 
(7). The Kramers-Moyal coefficients are the main ob- 
servables of a Markov process from which all the terms 
in the Kramers-Moyal operator will be determined. It is 
a well known theorem (Pawula theorem) of Markov pro- 
cesses that whenever the fourth order Kramers-Moyal co- 
efficient tends to zero all the the other terms with higher 
order derivatives tend to zero So the distinction 

between Markov processes which the Fokker-Planck de- 
scription when just the first two terms in the evolution 
operator in scale are important and the Markov processes 
which all the terms should be preserved is related with 
the coefficients too. Thanks to the detailed analysis car- 
ried over experimental data , the functional form of the 
first four Kramers-Moyal coefficients has been obtained. 
It is been observed that the fourth order conditional mo- 
ment tends to zero from which by invoking to Pawula's 
theorem the Fokker-Planck equation is canditated. In 
the mean time the present authors have shown that the 
Kramers-Moyal coefficients which are derived from the 
Yakhot modelling are extremely consistent with the ex- 
perimental observations It is interesting that the 
functional forms of the different coefficients (£)"([/, r)) 
are identically supporting the experimental observations 
0. The intermittency exponent of the structure func- 
tions can be derived from equation (7), C{p) — . It 
is easy to see that the ratio of different KM coeficients 
are controlled by the B parameter. As it is obvious when 
B — > oo K41 scaling is recovered and B ^ Q produces 
the extreme case of multiscaling related to Burgers inter- 
mittency 1^. The Kramers-Moyal's description of PDF 
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deformation in scale has been supported by exact compu- 
tation also for the compressible turbulence in high Mach 
number limit In that case the numerical values of 

the A and B terms is determined without any needs to 
numerical estimation. Reminding the original idea for 
suggesting the Markovian property of energy cascade in 
scale we take a step further and calculate the more gen- 
eral objects of the cascade, i.e. the unfused multiscale 
correlations. Assuming the Markovian nature of velocity 
increments in scale and the proposed form of evolution 
operator LKMiU,r), one in principle can calculate any 
correlation among velocity increments in different scales: 

Tn{x\ri,r2....,r„) - {U{n)U{r2)...U{rn)) = 
/ d(7(ri)...dC/(r„)C/(ri)...C/(r„)P(C/i,ri;;72,r2;...;f/„,r„) 

The joint probability P{Ui,ri\U2,r2] ■■■\UmT'n) can be 
calculated by advantage of Markovian property in terms 
of conditional probabilities, i.e., 

P(f7i,ri;C/2,r2;...;C/„,r„) = P(f/i, ri |C/2, r2) x 
xP(t/2, r2|C/3, r3)...P(t/„_i, r„_i|{/„, r„)P([/„, r„) (8) 

The conditional PDF of velocity increments can be writ- 
ten as a scalar-ordered operator as P(C/i, Ai |J72, A2) = 

( j^^ dXLKM(Uu\)) 

T(e_|_ ^ )^{Ui — 112)- So in the calculation 

of n-point multiscale correlation a series of conditional 
operators would emerge in the integrand of (^). When 
some of the coordinates coalesce the conditional oper- 
ator tends to a dirac delta function. The reduction of 
the conditional probability between the coalescing coor- 
dinates simplify the calculations. The only remaining 
conditional operator will be the probability of observing 
the typical velocity Ui increment between one subclass of 
fused points conditioned on observing the typical velocity 
increment U2 in the other subclass of fused points. Ex- 
plicitly we examine the behaviour of ^p4.g(Ai, A2) defined 
in 0, where Ai = ln(L/r) and A2 = ln(L/P). 

^P+,(Ai,A2) = (f/P(Ai)[/«(A2)> = j dUidU2S{Ui - U2) X 
X P (C/2,A2)(e-(^i-^^)^5cM(^i)[/f)(7| (9) 

We restrict the calculations to the GI invariant apprtox- 
imation so that neglecting the 0{urmsT' I L) operators in 
LiiAiiU, X). The crucial point in the above approxima- 
tion is that in the GI regime the Kramers-Moyal's coeffi- 
cients are scale independent so that all the scale depen- 
dence of conditional probabilities would reveal as a simple 
subtraction of the two logarithmic scales, i.e. Ai — A2 in 
the exponent. Because L^j^^(f/i)C/f = C(p)C^f , we will 
obtain the proposed form of the fusion rules in (2) with 
C(p) = (p+B) • other multiscale correlation function 
is also tractable under the same approximations. 
The Fusion rules first introduced through invoking 
two Kolmogorov type assumptions. The first one is the 



assumption of a scale invariant form for all the correla- 
tion functions in the inertial range. The second is called 
the universality meaning that when some arbitrary set of 
velocity differences in the correlation functions are fixed 
in the scale L , their precise choice will affect the cor- 
relation functions just as an overall factor. In terms of 
conditional averages the second proposition means that 

{U{nr\U{r2r) = S,{n)%,,{r2) (10) 

where it is assumed also that scale of r2 is in the order 
of integral scale while the ri is in the inertial range. The 
function <&p,g(r'2) isa homogenious function with a scaling 
exponent C„ — Cpi ^'^d is associated with the remaining 
n — p indices of T. Mathematically the above conditional 
correlation is easily verified 

(C/(ri)P|C/(r2)«) = Spin)UP/Sp{r2) 

In Yakhot modelling scaling hypothesis has been invoked 
from the very beginning of the theory when the rele- 
vant OPE terms has been chosen for closing the equa- 
tion governing over the generating function of longitu- 
dinal velocity increments. But we are showing that at 
least in the framework of Yakhot modelling the univer- 
sality proposition is the result of the Markovianity of the 
evolution of velocity increments in scale. On the other 
hand the necessary proof of Markovian property has been 
verified through the special scalar-ordered form of the 
conditional probabilities. This itself has arisen from the 
general invariances and scaling constraint of the Navier- 
Stokes Equation. So the universality condition in the 
language of multiscale correlation functions has in it's 
heart the very robust scaling invariance under the in- 
finite parameter scaling group ||]. We should empha- 
sise that non-universal effects of the large scale motions 
can also manifest through the scale dependent terms in 
the Kramers-Moyal operator. Still the general form of 
Universality assumption would be the leading behaviour 
while the 0{urmsr/L) will be the sub-leading corrections 
inducing the large scale effects [?] . 

Within the experimentally verified approximation that 
neglects the third and higher order KM coefficients |]ll|,0 , 
one can write the equivalent diffusion process in scale 
which dynamically gives the relation between velocity in- 
crements at two different scales. In fact approximating 
the KM equation with a Fokker-Planck evolution kernel 
can be interpreted as if velocity increment U is evoluted 
in "scale" A (logarithmic length scale) by the langevin 

equation 0, |^ = D(^^{U, X) + ^ D(^1{U, \)r){X) , where 
7](A) is a white noise and the diffusion term acts as a mul- 
tiplicative noise. Giving address to the Ito's prescription 
pO| the multipoint correlation function can be written in 
the form of a path integral as 

^(Ai,A2)= VUUP{Xi)m{X2)e-^>'i V^^H^TJ) x 

XP(C/2,A2) (11) 
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By a simple application of Bayesian rule probability den- 
sity in the outer scale A2 can also be written as a path in- 
tegral entering the information of the integral scale PDF 
which in a good approximation can be regarded as a 
Gaussian distribution Building up all the terms in 
a descriptive way the joint probability P{Ui, Ai; [/2, A2) 
is represented as a path integral over all possible paths 
between U{Xi) and C/(A2) transfering in an intermittent 
way all the information of the integral scale into the cal- 
culation. Without further attempt for calculating the 
multiscale correlation by path integral representation, we 
turn our attention to the Langevin dynamics instead. 
The resulting process is the well known Kubo [|lO| os- 
cillator multiplicative process. By using the Ito |10|] pre- 
scription one can deduce that 

Sx,U{x)^WiXuX2)6x,U{x) (12) 

The multiplier W{Xi , A2) can be easily derived in terms of 
ai and a2 and Wiener process at two logarithmic scales 
as: 

>V(Ai,A2) = e(""i(^i"^2)+v^[^(^i)"^(^2)l). 

Equation ( p^ ) encodes a simple cascade process. Cas- 
cade processes are simple and well known useful tools 
to describe the leading phenomenology of the intermit- 
tent energy transfer in the incrtial range. Both anoma- 
lous scaling exponents and viscous effects Q can be Un- 
der the framework of the derived KM equation we have 
shown that the equivalent cascade model from Navier- 
Stokes equation is derived by cutting the KM evolu- 
tion operator after third term. In fact we have given 
a recipe for finding the equivalent approximate cascade 
model from Navier-Stokes equation. Structure functions 
are described in terms of multiplier M^(Ai,A2) through 
'5'p('") = C'p([yV(r'/i)]P) where from the Langevin equa- 
tion pure power law arise in the high Reynolds regime 
{[W{r/L)]P) - {r/L)^^P\ In this approximation the scal- 
ing exponents would be — —pai + p{p — l)a2/2. 
From the direct calculation of Langevin equation one can 
easily find the behavior of the multiscale correlation func- 
tion Tp+q (r, R) . In the same framework, it is staightfor- 
ward to show that 

- {m^mimjW) - Spir)Sp+qiR)/SpiR) (13) 

The independence of multipliers in two diferent scales 
is allways assumed for the underlying cascade process 
otherwise the following relation would not be held. Our 
modelling equivalently encodes the following requirement 
by the obvious independency of increments in a Wiener 
process. Recently Benzi and coworkers Q analysed mul- 
tiscale correlation functions from high Reynolds exper- 
iments and synthetic signals. They have elegantly seek 



whether the fusion rules (|^) are compatible with the ran- 
dom cascade phenomenology. Their main result is that 
all multiscale correlation functions are well reproduced in 
their leading term ^ by a simple uncorrelated ran- 
dom cascade. In Yakhot modelling of the dynamics of the 
longitudinal velocity increments in scale all the above re- 
sults are recovered in the framework of Fokker-Planck 
approximation. It is also interesting also to seek the 
limiting behaviour of the Multiscaling correlation func- 
tion for Burgers turbulence which is tractable by taking 
the limit of _B — > in our formulation. Equation (3) 
shows that multiscaling correlation function will be in- 
dependent of the outer scale R which is consistent with 
our knowledge about Burgers turbulence jl^]. We think 
that preserving all the terms in KM equation would give 
the full information of cascade in length scale and this 
would give an answer to the question that wether there 
are other subleading processes acting for energy transfer 
from large to small scales. Preserving the GI breaking 
terms in the corresponding stochastic processes would 
also give the important unanswered question regarding 
the effect of uneven PDF's of velocity increments on the 
cascade models. 
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